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Abstract 

We report on the first results of self-consistent second order metric perturbations 
produced by a point particle moving in the Schwarzschild spacetime. The second 
order waveforms satisfy a wave equation with an effective source term build up from 
products of first order metric perturbations and its derivatives. We have explicitly 
regularized this source term at the particle location as well as at the horizon and at 
spatial infinity. 

1 Introduction 

The observation of gravitational waves opens a new window onto our universe and we also expect that 
the observation of gravitational waves will provides a direct experimental test of general relativity. 

The space mission LISA will primarily detect gravitational waves from inspiraling solar-mass compact 
objects captured by a supermassive black hole residing in the core of active galaxies. For these Extreme 
Mass Ratio Inspirals (EMRI) we use the black hole perturbation approach, where the compact object 
is approximated by a point particle orbiting a massive black hole. There are two nontrivial problems 
to consider in this approach: The self-force problem and the second order gravitational perturbations 
problem. Due to the self-force the orbit of the particle deviates from the background geodesic, i.e. the 
spacetime is perturbed by the particle itself. It is essential to take this deviation into account in order 
to predict the orbital evolution to the required order. For the headon configuration studied in this paper 
this was achieved in PQ. The gravitational self- force is, however, not easily obtainable for more general 
trajectories. 

We require the second perturbative order calculations to derive the precise gravitational waveforms 
to be used as templates for gravitational wave data analysis. In general, this computation has to be done 
by numerical integration. Hence, it is important to derive a well-behaved second order effective source. 
In this paper, we will focus on this later problem. 



2 Second order metric perturbations 

We consider second order metric perturbations (MP), g^ v = g^ v + h^} + hffi , with expansion parameter 
/i/M corresponding to the mass ratio of the holes and where g^ is the Schwarzschild metric. The 
Hilbert-Einstein tensor and the stress-energy tensor up to the second perturbative order is given by 

G,Ah v ] = G^[h^] + G$[hM) + G$[h (1 \hU] , T, v = TW + T^ + r#« , (1) 

where 

= — -h tiv>a ,a + h a ^.^y a - R a ^vh afi - -h;fj, u - -g^v{h\a aX - h-\ ,x ) , 
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G$[hW,hW] = R®[hV\hW]-±g^RW[hV\h<% 

ff (2)[l(l) 1,(1)1 1 1,(1) 1,(1) a/3 ,K(l)aP (h (l) , 1,(1) _ qlU) n 

-l(hW<*0 „ - lh^' a )(2h (1) ~h (1) ) + -h {1) h {1)a ^ ~-h {1) /i (1)/3;Q 

2^ i/ 3 2 «(w) a ^\a) + 2 ^ a -0 v 2 " 

On the other hand, the stress-energy tensor includes three parts, the first order stress-energy tensor 
which is the one of a point particle moving along a background geodesic; 

/ + °° J r M ^r" 

«W( a -*(r))H* (fc-, (2) 

where = {T(t), R(t),®(t),$(t)} for the particle orbit, the deviation from the geodesic by the self- 
force T^t' SpS> (See Ref. (Tj), and finally T^J h \ which is purely affected by the first order MP, 

where we have used the determinant g — g(l + h^) up to the first perturbative order. 



3 First order metric perturbations in the RW gauge 

Before considering the second order, it is necessary to discuss the first order MP, i.e., the Regge-Wheeler- 
Zcrilli formalism ^ [3J. The basic formalism has been given in Zcrilli's paper [3J, and it has been 
summarized in the time domain in [JJ[5]. In the following, equation numbers (Z:l), (Ll:l) and (L2:l) for 
instance, denote the equation (1) in [3J, [1] and [5], respectively. 

For the first order Hilbert-Einstein equation, G$ [h^] — 8ttT^J, we expand h^} and Tju) in ten 
tensor harmonics components. We then obtain the linearized field equations for each harmonic mode. 
For the even part, which has the even parity behavior, (— l) e , we have seven equations. We impose 
the Regge- Wheeler gauge conditions (RW), the vanishing of some coefficients of the first order MP: 
4 e ) = h[ e) = G = 0. 



We introduce the following wave-function for the even parity modes 
2r 



£(£+! 

where the suffix RW stands for the RW gauge. This function ?/>f™ n obeys the Zerilli equation, 

d 2 d 2 



(4) 



dt" d 



V" 



where r* — r + 2Mlog(r/2M — 1), the potential V[ vcn and the source <S*f™ n are given in Eqs. (Ll:l-2) and 
(L1:A.3). The reconstruction of the MP under the RW gauge have been expressed in Eqs. (L2:B.9-12). 

In the following, we consider a particle falling radially into a Schwarzschild black hole as the first 
order source. The equation of motion of the test particle is given by 



dR\ 2 = 2M\ 3 J_ + A _ 2AA 2 E _ / 2M\ dT(r) 



dt V R E 2 \ R J ' V R dr 



(6) 



where E and R are the energy and the location of the particle, respectively. The non-vanishing tensor 
harmonics coefficients of the energy- momentum tensor are A% m , and . Because of the symmetry 
of the problem we have only to consider even parity modes, i.e., described by the Zerilli equation. 

In the head on collision case, the MP in the RW gauge are C° (continuous across the particle). One 
can see this as follows. First, from Eqs. (Z:C7d) and (Z:C7e), we find that H2t m {= Hoi m ) and Ki m have 
the same differential behavior. Then, we note that d r Kt m ~ 9(r — R(t)), because the left hand side of 
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(U|) behaves as a step function near the particle. This means that Kg m (and also H^im) is C°. Using 
this, d r H 1£m ~ 0(r - R(t)) is derived from Eq. (Z:C7e), i.e., H Um is C°. (See Ref. [Sj.) 

Using the above fact, we can take up second derivatives of the function with respect to t and r. 

These quantities allow us to calculate the coefficients of the 5-terms in the second order source. 

Next, we consider the I = perturbations (1=1 modes can be completely eliminated in the center 
of mass coordinate system). In the RW gauge there are four non- vanishing coefficients of the MP, i?ooo, 
H\ oo, -^2 00 and -Koo- The gauge transformation has two extra degrees of freedom. We could choose the 
gauge so that -ffioo = -Koo = 0, i.e., the Zerilli gauge, but it is difficult to treat the second order source, 
since the MP are not C°. 

We instead consider here a new (singular) gauge transformation, chosen to make the metric perturba- 
tions C° and to obtain an appropriate second order source behavior. To derive this gauge transformation, 
we have also considered a regularization of the second order source at r = oo and the horizon at the same 
time 0. Note that we succeeded in choosing the gauge such that all of the above MP behave as C° at the 
location of the particle and vanish at r — oo and r = 2M. 



4 Second order Zerilli equation 

Since the first order MP contains only even parity modes, we can discuss the second order MP for the 
even parity modes only, i.e. in terms of the Zerilli function, 

X ™ {t ' r) = 2r + 3M (r 2 §i !C 2o(t,r)-(r-2M)Hi 2 o(t 7 r^j . (7) 

Here, we have considered the contribution from the £ = and 2 modes of the first order to the 1 — 2 mode 
of the second order since this gives the leading contribution to gravitational radiation. We also choose 
the RW gauge to second order. This Zerilli function satisfies the equation, Zf ven xf (i, r) = iSfo(*> r ) with 



87rV3(r - 2 M) 2 d 87r(r-2M) 2 d Sv^tt (r - 2M) d 

3(2r + 3M) d't W ^ r > + 2r + 3M Ft 2 ^"' 3 W 



4V2i7r(r-2M) 2 <9 (i) 8 V2 in (r - 2 M) (5 r - 3 M) M (1) 

I 20 V*; r ) ~ ~ , r . 2 -^20 \Pi r ) 



2r + 3M dr M v ' r(2r + 3M) 2 



8V3wr(r-2M) 2 d (0) 32^3^ (3M 2 + r 2 )(r-2M) g(0) 

3(2r + 3M) dr 20(hr)+ 3r(2r + 3M) 2 



The functions B20 etc. are derived from the second order quantities, G^[h^\ h^], T^J h ^ (and T^t' SF ') 
by the same tensor harmonics expansion as for the first order. 

The delta function ^^(a; - z(t)) in T$ h} includes 6 ( ^(n - Q(t)) = T,i m Y im(^)Y; m (Cl(T)). We 
have considered only the contribution from the £ — and 2 modes of the first order perturbations. 
Consistently, we use only the three components, /i^ 2 ^l2m(^)^m(^( T ))) ^(fc2)^ 0m (^)^0m(^( T )) an d 
/i{,i 0) F 2m (fi)y 2 * m (O(r)). 

We may wander if there is any (5 2 -term in the second order source. The answer is "No". This is 
because in the case of the head on collision, the MP under the RW gauge are C°: R$[hW,hW] includes 
second derivatives and we need one more derivative to construct iSf . (h^) 2 is C° x C° and its third 
derivative yields C° x 5' and 9 x 5 as the most singular terms. On the other hand, T^ 2 ;' 1 ^ includes only 
C° x <5 terms. 

From Eq. ([8]), we obtain the second order source as 

SUt,r) = ^(t.rJ + M^t,,), (9) 

where (2 ' 2) 5f and (0 ' 2) 5| are the contribution from (£ = 2) • (£ = 2) and (£ = 0) • (£ = 2), respectively. 

Note that while the above source term is locally well behaved near the particle location, some terms 
diverge as r — ► 00. So, we need to consider some regularization for the asymptotic behavior [7]. 



3 We could regularize and fix the gauge for the second order source. One could also proceed by first fixing the gauge. 
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In order to obtain a well behaved source for large values of r, we define a renormalized Zerilli function 

by 

x 2 z (^) = xio(^)-X2o g ' (2 ' 2) -X2o s ' W2) ; (10) 

rcg,(2,2) _ V 2 ( d „ RW ,, A „HW, 



-g.(o,2) _ 1 r (r - 2 M) / / M / r \\ ^ 9 ]m ■ 



^ = ~ 32^ (^V 1 ~ ~ ln \2m) ) F 20o(*> r)—H^ Q (t, r) 



+6 ll-5 7 jC(i I r)^(t ) r) + 4 I 1 + 3 - In j 1 H$&(t,r)-^H$&(t,r) 

+5 (l-4^)i«i,r)|^^ , 

where the suffix RZ means some gauge choice as discussed in Sec. [3J Finally, The best suited equation 
to solve numerically for X20 is then 

i 2 ° VCn x!o(^) - S^(t,r); (11) 



5 Discussion 



In this paper, we obtained the regularized second order effective source of the Zerilli equation in the case 
of a particle falling radially into a Schwarzschild black hole. Using this source, we are able to compute 
the second order contribution to gravitational radiation by numerical integration. 

(2 S F) 

To be fully second order consistent we have to include the term T^ v ' which is derived from the 
self- force on a particle. The self- force for a headon collision has been calculated in pQ, and in a circular 
orbit around a Schwarzschild black hole in [8], but have not been obtained in the general case yet. 

To prove that there is no S 2 term in the second order source, we have used the fact that the first order 
MP in the RW gauge is C°. In the general orbit cases (including circular orbits), the first order MP is 
not C° in the RW gauge, but it is C° in the Lorenz gauge [§]. This gauge choice favors the study of the 
second order perturbations for generic orbits. 
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